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' Abstract: 



A phenomenological description of single transverse spin effects in hadron-hadron 
inclusive processes is proposed, assuming a generalized factorization scheme and 
pQCD hard interactions. The transverse momentum, k±, of the quarks inside the 
hadrons and of the hadrons relatively to the fragmenting quark, is taken into ac- 
count in distribution and fragmentation functions, and leads to possible non zero 
single spin asymmetries. The role of k^ and spin dependent quark fragmentations 
- the so-called Collins effect - is investigated in details in p"p — > nX processes: it 
is shown how the experimental data could be described, obtaining an explicit ex- 
pression for the spin asymmetry of a polarized fragmenting quark, on which some 
comments are made. Predictions for other processes, possible further applications 
and experimental tests are discussed. 



1. Introduction 

Understanding the kinematics and the dynamics of the quarks inside the hadrons 
is a challenging and interesting problem. In simple quark-parton models the quarks 
inside a fast hadron are assumed to be collinear, i. e. their momentum is parallel to 
the momentum of the hadron. This assumption works well for unpolarized processes, 
like A+B — > C+X which are computed, at large energy and momentum transfer p?, 
by the usual convolution of distribution and fragmentation functions with a hard 
perturbative elementary interaction, according to QCD factorization theorem [I]]. 
However, in polarized reactions, the contribution of transverse intrinsic momenta 
often turns out to be crucial in understanding the experimental results. 

In Ref. 0, we proved how the large single spin asymmetries 

da^ — da^ da^ — da^ 

ANixF > PT) = d^TdTi = ^d^> (1) 

measured in p^p — > ttX by the E704 Group at Fermilab ||, can be accounted for 
by introducing k± and spin dependences in the distribution functions of quarks 
inside the initial polarized proton, da^^ stands for the differential cross section 
E n dPa p ' p~*^ x /cPp.^; xf = 2pi/\/s, where pi is the pion longitudinal momentum 
in the pp cm. frame, and ^/s is the total cm. energy. The proton spin is "up" or 
"down" with respect to the scattering plane: f (J,) is the direction parallel (antipar- 
allel) to p p T x p n . Other sizeable single spin asymmetries have been explained or 
predicted in Ref. [|J. 

The original suggestion that the intrinsic fc^ of the quarks in the distribution 
functions might give origin to single spin asymmetries was first made by Sivers 
such an effect is not forbidden by QCD time reversal invariance provided one 
takes into account soft initial state interactions among the colliding hadrons 0. 
A similar suggestion for possible origin of single spin asymmetries was later made 
by Collins |5[], concerning transverse momentum effects in the fragmentation of a 
polarized quark. A first simple application of this idea was given in ]7[. 

Qiu and Sterman || |9| have proven that a generalized factorization theorem holds 
in QCD with twist-3 distribution and/or fragmentation functions, which take into 
account initial or final state interactions. Recently ]l(J they have computed - using 
a simple model for a new twist-3 correlation function - single spin asymmetries for 
p^p — > ttX and p^p — > ttX processes. Their approach establishes a sound theoretical 
ground for simple phenomenological applications: a new quantity is introduced for 
valence u and d quarks which represents a correlation between these quarks inside 
a polarized proton and an external gluonic field; a simple parametrization of the 
new function is fixed by fitting data and some predictions for other observables are 
given. 

In this respect the computation of Ref. [TIJ is similar to that of Ref. 0: also 
in the latter a new quantity was introduced for valence u and d quarks, a quantity 
which requires some initial state interactions between the colliding hadrons in order 
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not to be forbidden by time reversal invariance, and the assumption that a factor- 
ization scheme still holds; a simple parametrization was given and fixed by fitting 
the experimental data, and predictions were subsequently made for other processes 



The main difference is that in Sivers and our approach 0, [|, [5| - based on 
the QCD improved parton model - the new quantity, called A N f a / p i (x a , k± a ) or 

21^ (x a , k± a ), has a partonic interpretation: it is the difference between the density 
numbers f a / p i (x a , k± a ) and f a / p i (x a , kj_ a ) of partons a, with all possible polarization, 
longitudinal momentum fraction x a and intrinsic transverse momentum k± a , inside 
a transversely polarized proton with spin f or |: 

A N f a / P i (x a , k Xa ) = f a/p i (x a , k ±a ) - f a/pl (x a , k ±a ) (2) 

fa/pi i^ai k±a) fa/pi {x a , fc_|_ a ) 

where the second line follows from the first one by rotational invariance. Notice that 
A N fa/pi {x a , k±a) vanishes when k± a — > 0; parity invariance also requires A N f to 
vanish when the proton transverse spin has no component perpendicular to k± a , so 
that 

A N fa/pi (x a , k ±a ) ~ k ±a sin a (3) 

where a is the angle between fcj_ a and the f direction. 

A N f by itself is then a leading twist distribution function, but its k± dependence, 
when convoluted with the elementary partonic cross-section, results in twist-3 con- 
tributions to single spin asymmetries. This same function (up to some factors) 



has also been introduced in Ref. |TT| - where it is denoted by ffp - as a leading 



twist T-odd distribution function. In Ref. [O the relation between f^p and twist-3 



correlation functions is explained and an evaluation of single spin asymmetries in 
Drell-Yan processes, originating from a gluonic background and the so-called gluonic 



poles, is given. The exact relation between A / and f^ is discussed in Ref. [JT3 



A function analogous to A N f a / p i (x a , k± a ) can be defined for the fragmentation 
process of a transversely polarized parton 0, giving the difference between the 
density numbers D h / a x(z,k±h) an d D h / a i (z, k± h ) of hadrons h, with longitudinal 
momentum fraction z and transverse momentum k±h inside a jet originated by the 
fragmentation of a transversely polarized parton with spin f or J,: 

A N D h/a] (z,k ±h ) = D h/a i(z,k ±h ) - D h/a i(z,k ±h ) (4) 
= D h/a i(z,k ±h ) - D h/a i(z,-k ±h ) . 



A closely related function is denoted by in Refs. []TTj, [14| and its correspondence 
with A N D is discussed in Ref. [|13[]. Again we expect 



A Dhjai {z, kj_h) ~ kj_h sin (3 (5) 
where (3 is the angle between k±h and the j direction. 
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We first discuss here the role of the above functions, (|2[) and (f|), in gener- 
ating single spin asymmetries in the large pt inclusive production of particles, 

+ B — > C + X, within a phenomenological QCD factorization model; we give 
explicit expressions for the single spin asymmetries taking into account the leading 
order contributions of parton transverse momentum (Section 2). In Section 3 we 
assume that only Collins effect is active - the analogous work for Sivers effect was 
performed in Refs. P, U - and show that it can explain existing data on p^p — ► nX, 
with some possible problems at the largest xp values; we obtain an explicit expres- 
sion for the function given in Eq. (|4|). Such function is then used to predict other 
single spin asymmetries in Section 4, while further comments are made and future 
applications, with attention to planned experiments, are discussed in Section 5. 

2. Single spin asymmetries in QCD parton model 

We write the differential cross-section for the hard scattering of a polarized 
hadron off an unpolarized hadron B, resulting in the inclusive production of 
a hadron C with energy Eq and three- momentum p c , A^ + B — > C + X, in a 
factorized form, as 



da 1 = 



E c da A]B ^ cx 



d 3 p, 



c 



= \y, Y, l^f 2 d 2 k ± ^ kli< ek ±0 ( 6) 

a,b,c,d X a ,X' a ,X b ,X c ,X' c ,X d ,X c 
Px a A X> a fa/Al(Xa, fc±«) A/b(z 6 , k ±b) M^.^ My c ,X d ,X' a ,X b ^A (*> k ^c) ■ 

If we choose the z-axis as the direction of motion of A^ and xz as the scattering 
plane, then the | direction is along the y-axis. 

The above expression is proven for collinear configurations, according to the QCD 
factorization theorem: we assume here its validity also when taking into account the 
parton transverse motion. There is no demonstration of the factorization theorem 
in such a case, and Eq. @ has to be considered as a plausible phenomenological 
model. 

In our expression (||) we have explicitely shown the k± dependence in the par- 
ton distribution and fragmentation. The usual number density of partons a inside 
hadron A is given by 



fa/A(Xa) = J d 2 k ±a f a/A l(x a ,kj 



Notice that while f a /A( x a) cannot depend on the parent hadron spin, f a /Ai( x a, k± a ) 
can, provided some soft initial state interactions are taken into account 0. ||. 

Px\i (x a ) is the helicity density matrix of parton a inside the polarized hadron 

A. The M\ c ,\ -,\ a) \ 's are the helicity amplitudes for the elementary process ah — > 
cd\ if one wishes to consider higher order (in a s ) contributions also elementary 
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processes involving more partons should be included. D x c ' A c , [z, k±c) is the product 
of fragmentation amplitudes for the c — > C + X process 

where the ¥) stands for a spin sum and phase space integration of the undetected 
Tx,\ x 

particles, considered as a system X. The usual unpolarized fragmentation function 
Dc/ C {z), i.e. the density number of hadrons C resulting from the fragmentation of 
an unpolarized parton c and carrying a fraction z of its momentum, is given by 

Dc/ C (z) = \ £ / d 2 k ±c I%X(z, k ±c ) . (8) 

A C ,A C 

We shall neglect in Eq. due to the limited Q 2 range of its application, the 
(unknown) Q 2 scale dependences in / and D; the variable z is related to x a and Xb 
by the usual imposition of energy momentum conservation in the elementary 2 — > 
2 process, which reads, in the collinear case, z = — (x a t + Xbu)/x a XbS, where s,t,u 
are the Mandelstam variables for the overall process AB — > CX, whereas s, i, w are 
those for the elementary process a& — > ccf. A similar expression holds when taking 
k± into account. The elementary process amplitudes are normalized so that 

do 1 1 ^-^ . * 2 

The cross-section do^ is readily obtained from Eq. @ by changing j into |; 
in absence of fej_, collinear configurations and helicity conservation in the pQCD 
elementary processes do not allow any single hadron spin dependence and it would 
result do^ = do^ (T^]. With non zero k± instead, spin dependences might still 



remain in f a /j\t( x aj k± a ) and/or in the fragmentation process of a polarized quark. 

We specialize now Eq. (§) to the case of p^p — > nX processes; we keep only the 
leading contributions in k±, i.e. we consider parton transverse motion only in those 
functions which would otherwise be zero [A N f a / p t (x a , 0j_ a ) = A N D h / a i(z, 0±h) = 0)] 
or integrate to zero, and neglect it elsewhere. From Eq. ffi), after some straightfor- 
ward but lengthy algebra, one then obtains, for p^p — > ttX: 



da r -do 1 = J 



X 



a,b,c,d ' 



1TZ 

„ ^£j.ab— >cd 

J d 2 k ± A N f a/ ^(x a ,k 1 _)fb/p{x b ) — -* — {x a ,x b , k±) D n/c (z) (10) 



dt 

+ J d 2 k' ± P a 'v ] f a/p (x a ) f b/p (x b ) A NN o ab ^ cd (x a , x b , k' ± ) A N D 7r/c (z, k' 
Eq. (|10D is worth some comments. 
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The second line of Eq. (|TU]) represents the Sivers effect: this expression of 
da^ — da^ is the one used in Refs. and || to fit the data [Q] on An, Eq. (TJ), 
and to obtain an explicit functional form for A N f u / p and A N fd/ p , Eq. (^). 

The third line of Eq. ( |i0| ) represents the Collins effect: this expression of 
dcr^ — da^ is the one we are going to use in the sequel to fit the data on A N 
and to obtain an explicit expression for some A N D 7T / C , defined in Eq. (^). 
Notice that 

A N D 7T/c (z,k' ± )=2lmD+J- c , (11) 

where [see Eq. ([?])] 

DtTc=t^ Vx**Vx xi - (12) 



j x,x x 

is a purely imaginary quantity, by parity invariance. 

pa/p T _ 2i p a J_ v _ is the polarization of parton a inside the transversely polarized 
proton pf; it is twice the average value of the f component of the parton spin. 
In this part only collinear configurations are taken into account. P a / pT might 
depend on x a , but we will assume that, at least for large x a valence quarks, it 
is constant. The product P a / pT f a / p {x a ) is the transversity distribution, often 
denoted by ATq(x a ) or hi(x a ) fll6| . 

Ajvjv<7(x a , Xb, k' ± ) is a double spin asymmetry for the the elementary process: 



A NN a(x a , x b , k' ± 



dt dt 

k E [M + ^M1^ -M+^M*-^] • (13) 



32 7T s , . 



The division of Eq. ( ^0|) into two separate pieces, corresponding respectively 
to spin and fc^ effects in distribution and fragmentation functions, is not due 
to keeping only leading terms in k±; the structure of Eq. fllCf ) - within its 
range of validity - remains unchanged also at higher order in k±, and one 
should only add the appropriate k^ dependences in all terms. 



3. Single spin asymmetries in p'p — > nX and Collins effect 

In Refs. and || we have considered A N f as the only possible source of single 
spin asymmetries in p^p — > irX processes and have obtained an explicit expression 
for it, by fitting the experimental data. This expression of A N f has been used in 
Ref. || to compute the values of An for p^p — > nX, in agreement with data, and 
to predict An for p^p — > KX, not yet measured. 

We consider here A N D as the only possible source of single spin asymmetries 
and see if we can get an equally good description of the data, obtaining an explicit 
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expression for it, to be used in other processes. We will find that this is possible, 
with some difficulties, and shall comment on it at the end of this Section. 

Of course, both effects - in distribution and fragmentation functions - might 
be at work at the same time and be in different proportions responsible for the 
observed asymmetries. By considering only one effect at a time we can obtain 
the maximum possible values for A / and A N D; such functions could then be 
used in processes where only one of them can be at work to produce single spin 
asymmetries, and see whether or not they give results in agreement with data. For 
example, in prompt photon production, A^B — > 7X [|J and Drell-Yan processes, 
A^B — > fi + fi~X, only A / can be active, whereas in semi- inclusive DIS, ip^ — > 
£ttX, only A N D should contribute, because initial state interactions in lepton-proton 



interactions are suppressed by powers of a em p7| . More data are expected in the 
future, from running or planned experiments at HERA, Jefferson Lab and RHIC; 
a richer and more precise amount of experimental information might also allow a 
simultaneous determination of both A / and A N D by using the full equation fllOl). 

There is another reason to explore here the effects of A N D alone. As we said, the 
function A N f is bound to be zero if the incoming hadrons are treated as free plane 
wave states |6|, |I~2| ; this conclusion may be avoided by invoking soft initial state 
interactions, different from those taken into account in the proof of the factorization 
theorem, together with the assumption that one can still use the factorized form (§) 
(see also Ref. fll8f ). The same does not apply to the fragmentation process, which 
requires anyway final state interactions; it might be - and this has to be determined 
experimentally - that the Collins effect is the main origin of single spin asymmetries. 

We consider then the last line only of Eq. (PUD to compute A N , Eq. ([[]), for 
p^p — > 7r ±,0 X processes: 

da^-dat = J^^d 2 k ± (14) 

a,b,c,d 

x P a/p1 fa/ P i (x a ) f b / P (x b ) A NN a ab - cd (x a , x b , kj_) A N D„ /c (z, k ± ) 



and 



da^ + da 1 = 2da unp = 2J2 I (15) 

a,b,c,d 

X fa/p(Xa) fb/p(Xb) —~ (X a , X b ) D n/c (z) . 

Let us consider Eq. flip) and notice that A N D W / C (z, k±), Eq. (f|), is an odd 
function of k±: this means that we cannot neglect the k± dependence of the other 
terms, and have to keep into account k± values in AjVAr<J a6 ^ cd . Eq. (0) can then 
be written as 

da^-da^ = £ J^^d 2 k ± P a ^ f a/p ,(x a )f b/p (x b )A N D n/c (z,k L ) 

a,b,c,d 

x [A NN &(x a ,x b , k ± ) - A NN a(x a ,x b ,-k±)] (16) 
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where now the integration on k± runs only over one half-plane of its components. 
The only unknown functions in Eq. (^) are P a /f T and A N D 7T / C (z, k±), which, at 
least in principle, are then measurable via the single spin asymmetry A^. 

In order to perform numerical calculations, we make some further assumptions, 
following Refs. and Q. Our first assumption is that the dominant effect is given 
by the valence quarks in the polarized protons. That is, we assume P a / pT to be non- 
zero only for valence u and d quarks: it is natural to assume a smaller polarization 
for sea quarks which, moreover, do not contribute much to the production of large 
xp pions. Secondly, we evaluate Eq. ( |T6"D by assuming that the main contribution 
comes from k± = k ^ where, as suggested by Eq. (|5]), k° ± lies in the overall scattering 
plane and its magnitude equals the average value of (k^) 1 ^ 2 = k\ (z). This average 
value does in general depend on z. The z dependence of the transverse momentum 
of charged pions inside jets was measured at LEP ]nj. We have made a fit to their 
data points. The result is shown in Fig. 1 and corresponds to the curve 



M 



0.61 z - 27 (1 - z) 



0.20 



(17) 



with M = 1 GeV/c 2 . 

The residual z dependence in A N D^/a not coming from k\_ is taken as a simple 
power behaviour, so that 



A N D(z,k° ± ) 



M 



N c z ac {l-zf c 



(18) 



where N c , a c and f3 c are free parameters. 

The k± integration then produces the simple expression: 



d 2 k ± A N D 7T/c (z, k_ 



A^ab^cd 
NN<? 



(k 



M 



N r z c 



NN a (k ± 



A~ab— *cd / 
NN& ( 



A±.ab^rcdi 
NNCT ( 



-k_ 

■kV 



(19) 



where (z) is given by Eq. (pTTf) . 

Exploiting isospin and charge conjugation invariance and assuming that a large 
z pion is mainly generated by the fragmentation of a quark which can be a valence 
quark for the pion, one can express all functions (|18D - one for each flavour c - in 
terms of a single one, 



A N D val (z } kl) = N 



M 



z a (i - z y 



(20) 



In fact 



A N D n+/u (z,k° ± ) = A N D 7r+/d -(z,k° ± ) 
A N D 7r - /u (z,k ( i) = A N D 7T - /d (z,k ± ) 



A N D val (z, k° 



(21) 



7 



and 

A N D 7T o /u (z,k° ± ) = A N D 7T o /d (z,k ± ) 
= A N D 7T o /u (z, fc° ) = A N D^ 0/S {z, k° ± ) = l - A N D val (z, fc° ) . (22) 

Notice that A N D w / c can be different from zero only for polarized quarks, i.e. 
quarks resulting from processes ab —>■ cd for which A^^a is not zero. These are 
the processes: qq — > qq, qq' — > qq', qq — > qq, qq — > qq and qg — > qg, and the 
corresponding Ann& can be computed from Eq. ([131). 

As initial polarized partons, like we said before, we consider only u and d va- 
lence quarks. The value of P g / pT is taken as a free parameter: given the overall 
normalization factor iV contained in A N D va i only the relative value 

P u/d = P u ^/P d ^ (23) 

is a meaningful free parameter. P u l^ is assumed to be 2/3. 

In the unpolarized cross-section we take into account all leading order QCD pro- 
cesses involving quarks and gluons. Expressions for the unpolarized partonic distri- 
butions and fragmentations can be found in the literature; we have chosen, as in Ref. 
the MRSG parametrization |2(J for the partonic distribution functions, and the 



parametrization (BKK1) of Ref. [|2T| for quark fragmentation functions into pions. 
These fragmentations functions are sufficiently well established and constrained by 
data; other sets available in the literature have similar large z behaviours and do 
not change significantly the quality of the fit and the conclusions which follow. 

Eqs. and ([15]) contain now known functions - the unpolarized distribution 
and fragmentation functions - computable elementary dynamics - da/dt and Aw^a 
- and unknown functions - A N D n / c and P a / pT - parametrized in a simple way, Eqs. 
(ID3)> (|20|)-(|22|) and (|23|), in terms of the four parameters N, a, (3 and P u /d- 

We have computed A N , Eq. (0), as a function of x p in terms of these parameters 
and have compared with the existing data for p^p —>■ n^' processes 0; in our 
computation the pt of the produced pion in the pp cm. frame has been fixed to 
Pt = 1.5 GeV/c, the average transverse momentum of the data. We obtain that the 
A N D va i function which best fits the data is given by Eqs. (|T7|) and (p0|) with: 



N = -0.22 

a = 2.33 for z < 0.97742 (24) 

P = 0.24 

and 

A N D va i(z,k° ± ) = -2D vai (z) = -2x1.102 z'^l-z) 1 * for z > 0.97742 , (25) 
where z = 0.97742 is the value at which \A N D va Az, k\ )| = 2D va Az). D va Az) is 



taken from Ref. f21fl . The best fit value of P u /d turns out to be 

P u/d = -0.76 . (26) 
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The corresponding fit is shown in Fig. 2. The fit is qualitatively good and satis- 
factory - one should not forget that we are considering an extreme situation, with 
only Collins effect taken into account and a simple k± dependence in A N D(z,k±); 
however, the last points at large xp values seem difficult to fit. In our computation 
A N — > when xp — > 1, due to the fact that k° ± (z) — > when z — > 1, Eq. ([17|) . We 
can obtain the best results only by letting A N D(z, ) to be as large as possible 
when z — > 1, i.e. we have to saturate the necessary inequality: 

\A N D(z,k° ± )\<2D(z). (27) 

In Fig. 3 we plot the ratio A N D va[ (z, k° ± )/2D va i; it shows the elementary left- 
right asymmetry, in the fragmentation of a transversely polarized valence quark, 



_ D„ /q i(z,k ± ) -D n/ql (z,k° ± ) 
D n/q r(z,k Q ± ) + D n/ql (z,k° ± ) ' 



which is needed to fit the data. 

I-A^/J has to reach the maximum value 1 at large z; smaller values would give 
too small values of A N at large xp. 

This was not the case when relying only on Sivers effect @, ||]. In Fig. 4 we 
plot the value of the left-right asymmetry in the splitting of a transversely polarized 
proton into a valence quark, 



j^N _ /<?/pT (X, fc°) f q /pl (X, ) _ A N f q / p (x, k° ± j 

./, ,(.'•■/.•") • ./;,„. ( --./.•°) " V; , 



as obtained from Refs. [§]] for u and d valence quarksQ. 

By comparing Figs. 3 and 4 we see that a small elementary left-right asymmetry 
Au,d/ P i n the proton distribution alone allows to fit the data on 0, U; instead, 
one needs a much bigger (at large z) left-right asymmetry A^ q in the quark frag- 
mentation alone in order to fit the same data. Moreover, the fit obtained here with 
the large A^, q might show some difficulty at the largest xp values. If such a dis- 
crepancy should be confirmed by further data it would be a significative indication 
towards the importance of spin and k± dependences in distribution functions. 

We also notice that in the simple model proposed by Collins for A N D 0, one 
finds that A^, q — > when z — > 1; such model could not explain the observed values 
of A N . 

The reason for the different results due solely to A N f [0, Q or A N D (this paper) 
can be qualitatively understood by looking again at Eq. ( |T0| ) and it turns out to 
be very interesting, as it is due to the different dynamical contributions in the two 



1 Notice that, due to a mistake in the writing of the paper, the values of N a given in Eq. (9) of 
Ref. are wrong and have to be multiplied by a factor 4; the fit to the data, Fig. 1 of Ref. Q], 
remains the same. 
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cases: let us compare Eq. ( |T9"D with its analogue for Sivers contribution, obtainable 
from the second line of Eq. (|T0"1) ||]: 



J d 2 k ± A N f a/p i(x a ,k 



da 



ab^cd 



M 



N n z aa 1-z 



Pa 



dt 
da ab ~* cd 



da 



ab^cd 



-1 



dt 



J "aft— >cd 

dt 



(30) 



In the expression for A at Eqs. ( |T9"D or (|30| ) have to be divided by the unpolarized 
cross section which contains /, D and da /dt. Now, it turns out that not only A NN a 
is smaller than da/dt, but also the cancellation between the two terms in squared 
brackets is stronger for A NN a, Eq. ([19]), than for da/dt, Eq. ([30]). This explains 
why fitting the data requires a much bigger value of \A N D\/2D than \A N f\/2f . 

We finally comment on the value found for P u /d which might be surprising: in 
SU(6), for example, one has P u ' p] = 2/3 and P d ' p] = -1/3, so that P u/d = -2. We 
find here that, indeed, the u and d valence quark polarizations inside a polarized 
proton have opposite signs, but the d polarization is, in magnitude, slightly bigger 
than the u polarization, Eq. (p6|). 



4. Application to other processes, p^p — > nX and p 



P 



KX 



We now apply the A N D vai (z, k° ± ) we have obtained by fitting the E704 exper- 
imental data on pion production, to predict single spin asymmetries in other pro- 
cesses. 

Recently the E704 Collaboration at Fermilab has presented results on single spin 
asymmetries for inclusive production of pions also in the collision of transversely 
polarized antiprotons off a proton target f22|. The kinematical conditions are the 
same as in the case of polarized protons ||. Within our model the connection 
between single spin asymmetries with polarized protons or antiprotons is very simple: 
since only valence quark contributions are taken into account, we just have to exploit 
the charge conjugation relations (|21|) and (|22"D for the A N D and similar ones for the 



pi/v j n particular, this means that we should expect 



A N (p i p -> vr ± ) ~ A N (p^p -> 7T T ) 



A N {fp -> tt°) ~ A N (p1p -> tt°) . (31) 



Our results are shown in Fig. 5, and indeed they agree with Eq. (|3l|) . The 
experimental data are in very good agreement with our predictions for 7r + and 
somewhat smaller than our predictions for 7r° and ir~ , although the overall agreement 
is good. The expectations ( |31~] ) are essentially based on charge conjugation invariance 
and the dominance of certain elementary dynamical contributions and it is difficult 
to see how they could not be true; we expect that further data both on p^p — > nX 
and p"p — > nX should lead to a better agreement with Eq. (pl|). 

We consider now the case where charged and neutral kaons are detected in the 
final state, p^p — > KX; the kinematical conditions are the same as for the other 
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processes. Lacking a determination of A N D(z, k° ± ) for the fragmentation of quarks 
into kaons we assume that the A N D val are the same for K and 7r production, up to 
a different normalization factor which we fix by imposing, for valence quarks, 

A N D K/q (z, fe° ) (Nk)_ 

where (N h ) = f$ dz D h / q (z). We have used unpolarized kaon fragmentation function 
sets from Refs. [f2T], |23|, |25] and have computed the appropriate normalization 



factor (Nk) /(N*) for each set. 

We should keep in mind that here only large x valence quarks from the initial 
polarized proton are taken into account: this means that contribution from strange 
quarks from the polarized proton sea are completely neglected and the production 
of kaons is through non-strange quarks only. The eventual role of strange polarized 
quarks (and antiquarks) could only be studied within a more refined model for the 
transversely polarized distribution functions, h\{x) = P q /^ f g / p (x), where contribu- 
tions from q = s and s are included. However, we have checked that, even assuming 
a large value of P s / pT and taking A n Dk/ s = A N D va i, the contribution of strange 
quarks to An is negligible at large xp, as f s / p (x a ) is small at large x a values. 

Fig. 6 shows the results we obtain for An in the production of either charged 
or neutral kaons, (K + + K~)/2 and (K° + K°)/2 respectively. The positive values 
of An refer to charged kaons, while the negative ones to neutral ones: the values 
obtained are large for all sets of unpolarized fragmentation functions, whose features 
are further discussed in Ref. H|. Notice that the different sets of fragmentation 
functions give qualitatively similar results, contrary to the findings of Refs. |3Q, 
where, for neutral kaons, the results show a strong dependence on the choice of the 
fragmentation functions. 

Let us stress once more that these results for K production are rather qualitative 
and based on the simple assumption about the form of A N Dk/ q (z, k ^), Eq. 
nevertheless, they clearly show how An might be large and measurable also for 
p^p — > KX processes and its measurement should give more information about the 
Collins effects in the fragmentation of a polarized quark into a kaon. 

5 Comments and conclusions 

Single spin asymmetries in inclusive hadron production can be explained within 
QCD; they are twist-3 effects which can be large and measurable in the pt range 
of few GeV/c. Several possible origins - leading to the introduction of new basic 
functions - have been discussed in the literature, but a clear and comprehensive 
understanding has not been reached yet: more theoretical and experimental work is 
still necessary. 

In this paper and previously in Refs. and Q we attempt a consistent phe- 
nomenological description and computation of single spin asymmetries based on 
QCD parton model and the introduction of fundamental spin asymmetries at the 
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quark level, both in the splitting of polarized protons into quarks, Eq. (g), and 
in the fragmentation of polarized quarks into hadrons, Eq. (^). The link between 
this approach and other approaches based on the introduction of parton correlation 
functions JED], |26| and/or gluonic poles [12[] has to be clarified. As suggested by 



Ratcliffe |26| the data might also be analysed using a general parametrization of the 
partonic correlators. 

In Ref. p| we showed how the so-called Sivers effect alone could account for data 
on An and we derived an explicit expression for the function A N f; here we have 
shown that also the Collins effect - with some more difficulties - can account for 
the same data, and we have obtained an explicit expression for the function A N D. 
We have also introduced a general formalism which includes both effects, Eq. (|T0|). 

Their relative importance has to be established by further experimental data; 
single spin asymmetries could be measured in the near future at Jefferson Lab, 



HERA [27| and RHIC. Of particular interest are those processes in which only one 
of the two effects is expected to be active: in semi-inclusive DIS, ip 1 * — > £hX, and 
in 7*p T — > hX initial state interactions are suppressed by higher powers of a em 
and single spin asymmetries should originate only from the fragmentation of a final 
polarized quark ||28|| . It would be very interesting to have data on such processes. 



On the other hand processes like p^N —>■ 7A and p*N — > fi + fi X could exhibit 
single spin asymmetries only due to spin and k± dependent distribution functions 



UK or quark-gluon correlations ||, [L^, Pgl ; also data on these processes would be of 
great help. 

A comprehensive phenomenological description and a physical understanding 
of subtle single spin effects, within perturbative QCD and a simple factorization 
model, seems indeed possible; more data should allow the determination of new 
basic distribution and fragmentation functions, allowing genuine predictions for new 
processes and a deeper knowledge of the nucleon structure and the hadronization 
mechanism. More theoretical work, towards a true generalization of the factorization 
theorem to include transverse momenta, is also needed and might be motivated by 
our phenomenological study. 
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Figure Captions 



Figure 1. The pion intrinsic average transverse momentum in a jet as a function of z. 
The diamonds are the data from Abreu et al. [19], the continuous line is our fit, as given 
by Eq. (17). 

Figure 2. Single spin asymmetry for pion production, p^ p — > it X. The data points are 
the E704 experimental single spin asymmetries [3] for ir + (diamonds), ir° (squares) and 
7r~ (triangles). The solid line is our best fit for 7r + , the dashed line for 7r° and the dotted 
line for 7r~, obtained under the assumption of Collins effect only. 

Figure 3. The ratio between lA^D^I and twice the valence unpolarized fragmentation 
function D, Eq. (28), as a function of z. Notice that \A N D\/2D = 1 for z > 0.97742. 

Figure 4. The ratio between \A N f va i\ and twice the valence unpolarized distribution 
function /, Eq. (29), as a function of x. The solid line refers to u quarks and the dashed 
line to d quarks. Notice that in both cases \A N f\/2f — ► for x — > 1. 

Figure 5. Single spin asymmetry for pion production mp^p — > it X. The data points are 
the E704 experimental single spin asymmetries [22] for ir + (diamonds), tt° (squares) and 
7r~ (triangles). The solid line is the result of our model (with Collins effect only) for 7r + , 
the dashed line for 7r° and the dotted line for 7r~. 

Figure 6. Single spin asymmetry, as predicted by our model (with Collins effects only), for 
the production of charged (upper curves) and neutral (lower curves) kaons, (K + + K~)/2 
and (K° + K°)/2 respectively. The solid lines correspond to the BKK1 set of unpolarized 
fragmentation functions [21], the dashed lines correspond to the BKK2 set [23], the dotted 
lines to the GR set [24], and the dashed-dotted lines to the IMR set [25]; see Ref. [4] and 
text for further details. 
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Figure 1: The pion intrinsic average transverse momentum kj_ in a jet as a function 
of z. The diamonds are the data from Abreu et al. [19], the continuous line is our 
fit, as given by equation (17). 
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Figure 2: Single spin asymmetry for pion production, p* p — > 7T X. The data points 
are the E704 experimental single spin asymmetry for 7r + (diamonds), 7r° (squares) 
and 7T~ (triangles) [3]. The solid line is our best fit for tt + , the dashed line for 7r° 
and the dotted line for tt~ , obtained under the assumption of Collins effect only. 
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Figure 3: The ratio between lA^D^I and twice the valence unpolarized fragmen- 
tation function D, Eq. (28), as a function of z. Notice that \A N D\/2D = 1 for 
z > 0.97742. 
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Figure 4: The ratio between \ A N f va i\ and twice the valence unpolarized distribution 
function /, Eq. (29), as a function of x. The solid line refers to u quarks and the 
dashed line to d quarks. Notice that in both cases |A Ar /|/2/ — > for x — > 1. 
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Figure 5: Single spin asymmetry for pion production in p' p — > n X. The data 
points are the E704 experimental single spin asymmetries [22] for ir + (diamonds), 7r° 
(squares) and (triangles). The solid line is the result of our model (with Collins 
effect only) for tt + , the dashed line for rr° and the dotted line for 7r~. 
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Figure 6: Single spin asymmetry, as predicted by our model (with Collins effects 
only), for the production of charged (upper curves) and neutral (lower curves) kaons, 
(K + + K~)/2 and (K° + K°)/2 respectively. The solid lines correspond to the BKK1 
set of unpolarized fragmentation functions [21], the dashed lines correspond to the 
BKK2 set [23], the dotted lines to the GR set [24], and the dashed-dotted lines to 
the IMR set [25] ; see Ref . [4] and text for further details. 
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